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Abstract 

It was shown by dos Santos et al. the importance of Alfven waves to 
explain the winds of Wolf-Rayet stars. We investigate here the possible 
importance of Alfven waves in the creation of inhomogeneities in the 
winds of early-type stars. The observed infrared emission (at the base 
of the wind) of early-type stars is often larger than expected. The 
clumping explains this characteristic in the wind, increasing the mean 
density and hence the emission measure, making possible to understand 
the observed infrared, as well as the observed enhancement in the blue 
wing of the H a line. In this study, we investigate the formation of 
these clumps a via thermal instability. The heat-loss function used, 
H(T,n), includes physical processes such as: emission of (continuous and 
line) recombination radiation; resonance line emission excited by electron 
collisions; thermal bremsstrahlung; Compton heating and cooling; and 
damping of Alfven waves. As a result of this heat-loss function we show 
the existence of two stable equilibrium regions. The stable equilibrium 
region at high temperature is the diffuse medium and at low temperature 
the clumps. Using this reasonable heat-loss function, we show that the 
two stable equilibrium regions can coexist over a narrow range of pressures 
describing the diffuse medium and the clumps. 

1 Introduction 

As demonstrated first by Lucy & Solomon (1970), the radiative momentum 
absorbed by UV spectral lines is able to initiate stellar winds, since the radiative 
line acceleration exceeds the gradient by a large factor. The first model to derive 
mass-loss rates (M) and flow speeds in good agreement with observations was 
that of Castor, Abbott & Klein (1975) (CAK). One of the major difficulties 
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presented by the radiation driven wind theory is the momentum problem in 
WR stars, which can be described using the ratio r\ = (M« 00 )/(L»/c), where 
Woo is the terminal velocity and L* the star luminosity. Barlow et al. (1981) 
found that in WR stars 77 ranges from 3 to 30. This means that there is about 
an order of magnitude more momentum in the wind than in the radiation field. 
It was assumed that every stellar photon transfers its momentum, hv/c, only 
once (single scattering) , but even with multiple scattering of the photons one 
obtains Mv^ > bL^/c. To get around the momentum problem, one cannot 
simply appeal to a larger luminosity, because the values that are used cannot 
be near the Eddington limit (Cassinelli & van der Hucht 1987). 

Following the suggestion that there may be appreciable magnetic fields in 
WR stars larger than 1000 G (Mahesvaran & Cassinelli 1988; Poe et al. 1989), 
it was suggested that the wind in a WRN5 star, for instance, can be driven 
by Alfven waves (see Hartmann & Cassinelli 1981). They assumed B=20,000 
G and a mechanical flux of Alfven waves of <& w = 1.1 x 10 14 erg/cm 2 s (this 
work did not take into account the contribution of the radiation pressure on the 
lines). 

As implied by the work of Willis (1991), an additional mechanism to 
radiation pressure may be required to initiate the high WR mass-loss, although 
thereafter the winds may be radiatively accelerated. In this context, dos Santos 
et al. (1993a, b) proposed a model for mass-loss in WR stars, where both a flux 
of Alfven waves and radiation pressure are considered. The model is a fusion 
of the Alfven wave wind model of Jatenco-Pereira & Opher (1989a, b) and the 
radiation pressure CAK model. In the model an effective escape velocity is used, 
which takes into account the CAK power index expressing the effect of all lines, 
possible nonsolar abundances, and the finite size of the star disk. Their work 
indicates that Alfven waves, acting jointly with radiation pressure, provide the 
necessary energy and momentum for the wind, with reasonable Alfven fluxes 
and magnetic fields. 

Early-type stars show superionization lines O VI, N V, H a and X-rays, that 
cannot be explained by the high temperature star. Applying the coronal zone 
model to the winds of early- type stars, Cassinelli & Olson (1979) derived the 
ionization conditions expected in the wind of Q Pup. The results of this study 
explain very well the persistence to low effective temperatures of the strong lines 
of O VI, N V, C IV and Si IV. 

Since Abbott et al. (1984), one knows that the observed IR emission is often 
larger than expected from a homogeneous wind. From that time it was pointed 
out that the clumping in the wind, increases the mean density and hence the 
emission measure. Clumping can also explain an observed enhancement in the 
blue wing of the H a line. The narrow absorption components are likely to be 
direct manifestations of dense clumps. Now, the existence of these clumps is 
largely known in many individual hot stars, and clumping may be important 
in all hot stars with winds (Hillier 1991; Robert 1994; Moffat & Robert 1994; 
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Massa et al. 1995; Brown et al. 1995; Moffat 1994,1996; Eversberg et al. 1996). 

In principle, it can be said that the series of papers by Owocki, Rybicki and 
Castor (Owocki k Rybicki 1984, 1985, 1986, 1991; Owocki et al. 1988; Rybicki 
& Owocki 1990), which contains numerical hydrodynamic calculations, show 
that the radiation driven winds are violently unstable and that the consequent 
shocks can explain the X-ray emission of early-type stars, moreover the clumping 
explains the infrared emission excess and the formation of narrow absorption 
components. These models qualitatively explain the hot stars wind structure, 
and the role of Alfven wave damping in these winds is an open question, as 
noted by these authors. In this sense, our model is an attempt to explain the 
inhomogeneities of hot star winds using thermal instability in the presence of 
Alfvenic heating. 

The propagation and transmission of magnetohydrodynamic waves through 
stellar atmospheres and winds has attracted considerable interest, because of 
its relevance to the questions of chromospheric and coronal heating and wind 
acceleration (Leer et al. 1982). In order to accelerate a wind efficiently, the 
waves must be able to propagate without much reflection or attenuation up to 
the sonic point, because any addition of momentum below that point essentially 
goes to increase the mass flux but not the asymptotic wind speed (Leer et al. 
1982). On the other hand, as commented by Velli (1993), waves reflected and/or 
having a nonlinear evolution in the lower atmospheric layers, contribute to the 
nonradiative heating through turbulent decay. 

We study here a mechanism to form condensations in the base of early- 
type star winds. The basis for this approach is the work of dos Santos et al. 
(1993a, b) - a wind acceleration model for Wolf-Rayet stars where Alfven waves 
act jointly with the radiation pressure. Considering the above scenario, with 
special attention to the clumping features, our intent is to study the relevance 
of a flux of Alfven waves in the hot medium present in the base of WR winds, 
in order to understand the formation of the clumping features. As the magnetic 
field is more effective in the base of the wind (see dos Santos et al. 1993a), we 
have there a wave flux, resulting in a model that could explain clumps in this 
region. We consider a thin hot corona atmosphere, which corresponds to the 
X-rays observed, with a temperature of <~ 10 7 K, and density of ~ 2 x 10 13 cm~ 3 
(see van der Hucht 1992). In principle, the thermal instability process presents, 
as a result, the clumps observed, and we call these clumps the cool atmosphere 
(~ 10 4 K and ~ 10 times denser than the hot medium) (for instance, Brown 
et al. 1995). Our final goal is to demonstrate the stability of the base of the 
wind that has "cool" (~ 10 4 K) clouds and "hot" (~ 10 7 K) intercloud medium 
coexisting at the same pressure. 



3 



2 Physical Mechanisms 

In general, astronomical objects are formed by self-gravitation. However, some 
objects cannot be explained by this process. For these objects the gravitational 
energy is smaller than the internal energy. In these cases, it is assumed that the 
internal pressure is balanced by the pressure of the external medium. These 
objects (that cannot be explained by self-gravitation) are formed from the 
medium by some kind of condensation process not involving gravitation. Parker 
(1953) argued that, if the thermal equilibrium of the medium is a balance 
between energy gains and radiative losses, instability results if, near equilibrium, 
the losses increase with decreasing temperature. Then, a cooler-than-average 
region cools more effectively than its surroundings, and its temperature rapidly 
drops below the initial equilibrium value. 

Following Lucy & Solomon (1970), for a given ionization potential, \r = hv r , 
the photoionization rate, T r is a function only of the radiation temperature, T r ; 
the collisional rate, r c , on the other hand, is determined by the electron density, 
n e , and the electron temperature, T e . The ratio of the two rates is 

IV 

r c 

here ay is the cross-section for collisional ionizations, v e the electron 
velocity, and a„ the photoionization cross-section. B6hm (1960) has given 
approximations for these quantities from which one derives 

r r 1/2 v 3 
r c n e 

where T r — T e , hv r /kT r > 1, approximately, and the units of x are electron 
volts. 

Applying these results to the ionization of C III ions, for instance, \ r = 47.9 
eV with logT e // = 4.5. Taking T e (in the hot atmosphere) ~ 10 3 T e //, and 
n e = 2x 10 13 cm~ 3 , one obtains jf- = 3 x 10 5 (for lower density the ratio is even 
higher), so that collisional ionization may be completely neglected. 

Considering a thermal instability in an isobaric regime (e.g., Field 1965) 
(internal pressure balanced by the external pressure), we looked for a set of 
physical parameters that, at equilibrium [H(T, n) — 0], show three equilibrium 
regions: one stable region representing the diffuse medium; one unstable 
region; and another stable region representing the condensations. The energy 
gains considered are: heating by photoionization-recombination, H r ; Compton 
heating, H c ; and Alfven wave heating, Ha- These gains are balanced by the 
following radiative loss processes: cooling via thermal bremsstrahlung, Hf,; 
inverse Compton cooling (this term is computed jointly with H c ); and collisional 
excitation followed by resonance line emission, H T \. 
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2.1 Bremsstrahlung Losses 

The total amount of energy radiated in free-free transitions, per cm 3 per sec, in 
the case of a Maxwellian distribution of velocities, is 

/27rfcT\ 1/2 2 5 7re 6 , 



H b = -1.42 x 1CT 27 Z 2 n z n e T 1/2 g ff 



H b = -A fc T 1 /2 n 2 . (i) 

The quantity gff appearing above is a correction factor required for precise 
results. Its value is generally about unity (Spitzer 1978), and Ah = 2.4 x 10~ 27 . 
Hereafter n e = n, the number density. 

2.2 Resonance Line Emission 

Raymond et al. (1976) calculated a radiative cooling coefficient for a low density 
gas, optically thin, with cosmic abundances, between the temperatures of 10 4 K 
and 10 8 K, which we adopt in this work. A good fit to radiative losses, in 
this temperature range, due to electron excitation of resonance transitions in 
common metal ions (erg/cm 3 s) is (Raymond et al. 1976; Mathews & Doane 
1990), 

Hri = a ~^r n 2 , (2) 

1 + bTi w 

with a = 1.53 x 1CT 27 , b = 1.25 x 1(T 9 , p = 1.2, and q = 1.85. 

2.3 Photoionization-Recombination Heating 

An approximate equation to express the residual heating due to radiative 
ionization followed by recombination, in erg/cm 3 s, is 

H r = a B (T)max [0, (hv)i - hv a - f3k B T/2] n 2 , (3) 
where {hv) il the mean energy of ionizing photons, is 

J2 n 3 hv 3 _ n H kT H + n L kT L 



{hv)i 



n i n H + n L 



til and uh are the numerical densities of photons of the low and high 
temperature regions, respectively, given by 
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oT 4 

riL = , J? with T L = T* , 

and 

[H r i(T H ,n H ) + H h (T H ,n H )]tR* 
n H = 



kT H c 

where a B = 2.60 x 10" 13 (T/10 4 ) 
and /i^ the ionization potential of hydrogen. In the above equations uh and 
Th are the density and temperature of a high temperature region, nj, and Tl 
are the same for a low temperature region, and tR* is the region thickness. Each 
recombination results in a loss of energy /3/csT/2 from the thermal energy of 
the plasma, with / w 0.43 (Mathews & Doane 1990). 

The recombination expression, eq.(3), is an approximated one. Equation (3) 
states that we have recombination only when the average energy of the photons, 
(hv)i, is sufficiently high such that ionization can occur, that is, when (hv)i is 
greater than hv§ + f3ksT/2 (i.e., the sum of the ionization potential plus the 
average energy of the electron that is liberated) . 

2.4 Compton Heating-Cooling 

We have to estimate the number and frequency of the photons acting in 
the immediate neighborhood of the star surface and anywhere in the star 
atmosphere. Taking into account the interaction between thermal electrons and 
the radiation field, photons with lower frequency, come from a cooler optically 
thick region, the stellar continuum. Their flux is crT* 4 , and then, for these 
photons, we have 

LbL rp 4 

- — = crJ* 
AttR l 2 

On the other hand, a hot region of thickness tR* (optically thin), at temperature 
Th causes heating in the medium via, principally, thermal bremsstrahlung and 
resonance line emission [Hb(T 7 n) and H r i(T,n)]. Hence, for these photons, 

LbH 



4ttR 2 h 



[H b {T H , n H ) + H rl (T H , n H )] tR* 



The complete expression for Compton heating and cooling is then 

LbL 



ti r = 



(4) 



which is similar to the expression usually adopted, for instance by Mathews 
& Doane (1990). In the above equations ks is the Boltzmann constant, ctt 
the Thompson cross section, a the Stefan-Boltzmann constant, n the number 
density, m the electron mass, T the temperature, Lb the bolometric luminosity, 
c the light speed, T* the stellar temperature and i?* the stellar radius. 
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3 Damping and heating from Alfven waves 



Alfven waves in a early-type star, whose winds are primarily radiatively 
driven, are subject to damping (or amplification) as described, for example, 
by MacGregor (1996). In this case, the dispersion relation for Alfven waves in 
a radiatively driven wind, is k 2 v A = oj 2 — iujiv (instead of k 2 v\ = ui 2 ), where 



I*(v) is the intensity of the photospheric radiation field, ul is a line rest 
frequency and kl is the line mass absorption coefficient. If dl*{v)/dv > 
then the Alfven wave is amplified, while if dl*(v)/dv < the Alfven wave 
is damped. Although, as MacGregor (1996) noted, "the presence of such 
radiatively modified Alfven waves in the flow has yet to be explored" . We apply 
in the present investigation, the dampings described below with their heatings. 

The damping mechanisms that we assume here were used before in many 
astrophysical objects: protostellar, late-type stars and solar winds (Jatcnco- 
Pereira & Ophcr 1989a, b); galactic and extragalactic jets (Opher & Pereira 
1986; Goncalves et al. 1993b); early-type stars (dos Santos et al. 1993a, b); 
broad line regions of quasars (Goncalves et al. 1993a, 1996); cooling flows of 
galaxy clusters (Friaca et al. 1997) and others. 

3.1 Nonlinear damping 

Parallel Alfven waves are purely transverse and there is no important linear 
damping. The damping that does occur is not linear and it arises from a beat 
wave (two circularly polarized parallel propagating waves) which contains a 
longitudinal field component and a longitudinal gradient in the magnetic field. 
This results in a nonlinear damping of both electrostatic and magnetostatic 
components i.e. transient time damping. 

Volk & Cesarsky (1982) derived an equation that represents the unsaturated 
Landau damping, in the case of nonlinear two-wave interaction, that can be 
written as 



where 3 is the energy density in waves normalized to the ambient magnetic 
energy density, B^/Stt (Lagage & Cesarsky 1983). Using 3 = s^sl anc ^ ^ = 7u ' 
we obtain: 



too = 



2irn L i' L dl*(v) 
3c 2 dv 






where £ = 5 — 10 and v s is the sound velocity. 
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3.2 Turbulent damping 

There are strong evidences favoring anisotropic, supersonic and compressible 
turbulence in WR winds. Since all WR stars observed intensively so far do 
behave similarly, and WR stars are extreme manifestations of winds in hot 
luminous stars, it is possible or even likely that all hot-star winds show the same 
basic phenomenon (Moffat et al. 1994 and references therein). A necessary (but 
not sufficient) condition that one is dealing with turbulence is that the Reynolds 
number be Re >> 1. For an expanding wind, with v w (r) the expansion speed 
and r the distance from the star, one has: 

rv w (r) r v w (r) 
He w w , 

^thermal <>mfp V s 

where v is the viscosity, l m fp the mean free path of the average particle in the 
medium. Thus, with typical WR wind values, where the observed lines form, 
Re is much higher than 1, so turbulence is likely to exist if there is a driving 
force. 

Hollweg (1986) considered a new hypothesis for the nonlinear wave 
dissipation of Alfven waves. The hypothesis is that the wave dissipates via 
turbulent cascade, or, this hypothesis concerns the consequences of the Alfven 
wave dissipation in terms of wave-particle interactions, where the required 
power at high frequencies is presumably supplied via turbulent cascade. Then, 
exploiting the similarity of Pb oc fc~ 5 / 3 and Kolmogorov turbulence in ordinary 
fluids, the plasma volumetric heating rate associated with the cascade is given 

ft-^, (6) 

where p is the mass density, (Sv 2 ) is the velocity variance associated with 
the wave field, and L corl is a measure of the transverse correlation length. 
A subhypothesis is that the correlation length scales as the distance between 
magnetic field lines, 

In spite of L co „ being a free parameter, the model comes close to the notion 
of a Kolmogorov-like cascade to small scales. The waves themselves are here 
regarded as the source of the heating. In this case (Sv 2 ) only includes the power 
associated with Alfvenic fluctuations. Similarly, L corr concerns the correlation 
length of the Alfvenic fluctuations. Finally, in terms of damping length, we have 

L t = L COII v A (Sv 2 ) , (7) 

(Hollweg 1986, 1987). 
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3.3 Alfven wave heating 

Data of the last 10 years show us that early-type stars can be separated in 
two groups: magnetic stars, with surface strengths of a dipolc or quadrupole 
magnetic field of B n(10 2 — 10 3 )G, n = 2, 3, ...7; and normal stars, with 
B w 1 — 100 G. The magnetic field strength increases towards the center of 
the star and in the core is « (0.1 — 10) x 10 6 G, depending on the stellar mass 
(Dudorov 1994; Bohlcndcr 1994). The origin of these fields is an open question, 
and two theories compete to explain it: dynamo and fossil theories (Moss 1994). 

Consider now a collapsing cloud. For the collapsing cloud the cross sectional 
area perpendicular to a magnetic field, A, is oc p~3 and B oc A^ 1 oc ; where p 
is the mass density of the gas and B is the magnetic field. The damping length 
in each case is L = iu/T (i.e., the ratio between the Alfven velocity and the 
damping rate). Knowing that va = B/y/4np, p(Sv 2 ) oc $> w /va, where <£% is 
the wave flux and Q w oc p 2 ^ 3 , we write the nonlinear Alfvenic heating as: 

Hru = * W /L nl = ^T nl oc %^ oc „-V2 ; (8) 
va v\ B l 

and for the turbulent Alfvenic heating, 

1/2 

H t = * w /L t OC ^BV2vM_ K 7/12 . (9) 

VA 

following equations (5) and (7). 

The sum of the contributions from Compton and inverse Compton, 
photoionization-recombination, bremsstrahlung and resonance line emission, is 
about 10~ 22 n 2 . We are adopting ur (the density of the hot atmosphere)= 
2 x 10 13 cm~ 3 , resulting 10~ 22 n 2 w 4 x 10 4 erg/cm 3 s. We then normalize the 
Alfvenic heatings using F nh F t = 10 2 - 10 5 erg/cm 3 s. So, 



and 



Hnl ~ Fnl [ 2 ^ 0x 1Q13 ) 



3.4 The overall heating/cooling behavior 

In order to make clear the relevance of each heating/cooling process in 
the overall balance we plot, in Figure 1, the module of the heating or 
cooling due to: resonance line emission; photoionization-recombination; thermal 
bremsstrahlung; Alfvenic turbulent and nonlinear heatings; and Compton 
interactions, as a function of temperature, in units of erg/cm 3 s. 
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The first characteristic we note from Figure 1 is the fact that resonance 
line emission is the most important cooling. In fact it dominates over all the 
other processes in this range of temperature. Another aspect is that we are 
using only the contribution of the photoionization-recombination processes that 
produces heating. This mechanism is not considered at temperatures higher 
than logT ~ 6.3. At these high temperatures it appears as cooling (see also eq. 
(3))- 

From the physics of Alfven wave heating, it is clearly not temperature 
dependent (eqs. (10) and (11)), beyond a dependence on the density which 
implicitly scales inversely to temperature, keeping P fixed. Then, as Alfven 
heating is proportional to n a (a = —1/2,7/12), at a given pressure, we have 
this heating proportional to T~ a , as can be seen from Figure 1. 



4 Results 

The complete heating-cooling function, H(T,ri), including the physical 
processes discussed above, is: 

H(T,n) = -\ b T*n 2 - n 2 + - f x 

1 + bTi rac z 

{ (T H - T) [ H b (T H , n H ) + H rl (T H , n H )] tR. + (T L - T) nT^) + 

f3k B T] 



as (T) max 



0, (hu)i - hv a 



n 2 + H A , (12) 



with Ha assuming the form of H n i and Ht, given by equations (10) and (11). 
All the constants in (12) are in cgs units. 

Figures 2 to 5 show the balance between energetic gains and radiative losses 
(from equation (12)), i.e., the equilibrium of H(T,n)(~ 0) in a logP — logT 
diagram. For this calculation we assume hr = 2.0 x 10 13 cm~ 3 , Tjj = 1.0 x 10 7 K, 
T L = T*= 3.0 x 10 4 K; R* = 12i? Q and t = 0.2. 

As we are forming clouds, via thermal instability, from the hot atmosphere 
(n H = 2.0 x 10 13 cm~ 3 and T H = 1.0 x 10 7 K), we performed calculations in 
order to find, for each temperature, the density that corresponds to the balance 
(H(n,T) w 0). From the isobaric instability criterion, we need the clouds and 
the hot medium coexisting at the same pressure. In the hot atmosphere the 
pressure is Pr = 2nnkTH ~ 5.5 x 10 4 dyn/cm 2 s. We then want to verify the 
possibility of forming clouds (~ 10 K, cool and dense atmosphere) at pressures 
near the characteristic pressure of the hot medium. 

In Figures 2 - 5 we assume a single value of Fa for the two Alfvenic heating 
terms, nonlinear (F n i) and turbulent (F t ). In Figure 2, we have the equilibrium 
diagram in terms of pressure and temperature adopting Fa = 10 2 . This figure 
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shows us the coexistence of the two-phase equilibrium at pressures lower than 
the characteristic pressure of the hot atmosphere (logPfj = 4.74). For this 
value of Fa, the nonlinear heating (in the case in which the Alfvenic heating 
contribution to the overall heating-cooling function is the nonlinear one) , as well 
as the turbulent heating, are insufficient to permit the coexistence of the hot and 
cool atmospheres, at appropriate pressures. Analyzing Figure 3 (Fa = 5 x 10 3 ), 
we observe that both cases (nonlinear and turbulent) are satisfactory, in order 
to permit cloud formation via thermal instability, since they reach the pressure 
desired. Finally, in Figure 4 (F A = 10 4 ) and Figure 5 (F A = 10 5 ), the cool and 
hot stable solution can be found at the pressure desired (the reference pressure 
and higher pressures) in the two cases. In addition, it is clear that turbulent 
Alfvenic heating does a better job than nonlinear heating, since the pressure 
range of the coexistence of both stable equilibria is bigger (see bottom panel in 
Figures 4 and 5) in the case of this Alfvenic heating. 

5 Discussion and Conclusions 

Our results can be discussed in terms of the efficacy of the Alfvenic heatings in 
forming condensations near the surface of WR stars due to thermal instability. 
This efficiency is included in the scaling factor Fa, that is equal F n i or F t in 
Figures 2-5. 

We consider in this study a heuristic derivation of the expressions for Alfvenic 
heatings. The nonlinear and turbulent Alfvenic heatings represent extreme 
opposite dependencies of these heatings on density. Comparing this behavior 
we have: turbulent heating (H t oc n 7 / 12 ) which deposits more energy when 
the density is higher; on the other hand, nonlinear heating (H n i oc nT 1 !" 1 ) 
deposits more energy when the density is lower. These behaviors also can be 
seen from Figure 1 in which cooler regions are denser than hotter ones (that 
figure was plotted for a fixed pressure (Ph), then, the nonlinear heating curve is 
a decreasing function of density, opposite to the case of turbulent heating). Due 
to the completely different behavior of these heatings, the results from each 
one are very different (see plots 2 to 5 with the equilibrium solution for our 
models). Despite the fact that the Alfvenic heatings do not work in the cool 
solution, as well as in the hot one, results with nonlinear heating produce the 
cool condensations at about 2.2 x 10 4 K, while the results with turbulent heating 
show cooler low temperature solutions (<~ 1.0 x 10 4 K). Noting also the way that 
each process scales with P, at fixed T, one can understand why the stable hot 
solution has a narrower range in pressure than the stable cool solution, in all 
figures of equilibrium (figures 2-5). 

In the work of dos Santos et al. (1993a), the principal emphasis was to 
determine the terminal velocity of the Wolf-Rayet star winds, using a model 
which had radiation pressure and Alfven waves driving the wind. The initial 
Alfven wave flux, required was w 5.6 x 10 12 erg/cm 2 s. Using this value for 
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the wave flux, we can estimate the damping length for the Alfven waves. As 
in the derivation of Alfvenic heatings (subsection 3.3), & w is equal to HaLa- 
Adopting, for example, the maximum pressure in which the stable two-phase 
equilibrium exists, in Figure 5b, Ph - 5x 10 5 dyn/cm 2 , the turbulent Alfvenic 
heating for the formed clouds is ~ 1.9 x 10 7 erg/cm 3 s. Then, 

5.6 xlO 12 6 

Lfi rnir .\ = ~ =- ~ 3 x 10 cm. 

t(mm > H t 1.9 x 10 7 

Taking now the minimum value for the pressure in the cloud in Figure 5b, 
P ~ 3.9 x 10 4 dyn/cm 2 , the turbulent Alfvenic heating is <~ 4 x 10 6 erg/cm 3 s. 
We then have 

<f> w 5.6 x 10 12 „ » 

L t (max) = TT — ~. TTTTT" — I- 4 x 10 Cm. 

{ > H t 4 x 10 6 

These values for the Alfvenic damping lengths can be understood as limits to the 
size of the formed clouds, since the cloud diameters (d c ) must be smaller than 
the damping lengths (3 x 10 5 < d c < 1.4 x 10 6 cm), in order to have turbulent 
damping of the Alfven waves effective in this cloud formation process. 

The Alfven flux adopted here was the necessary value in order to accelerate 
the wind to the observed velocity, and obtain the necessary momentum in the 
wind, with minimum magnetic field. The ratio between this flux and the total 
one, at the star surface is 

£ = ® w 
LjAuR 2 

Assuming typical values for (10 5 ' 5 L©) and i?» (12 Rq), we obtain £ w 0.04, 
which means that Q w is only a few percent of the total stellar flux at the star 
surface. (This flux of waves is lower than others used in the literature, for 
instance, Hartmann & Cassinelli 1981.) 

It is interesting to estimate also the magnetic field before (in the diffuse 
medium, B a ) and after (in the clouds, B c ) the condensation process. During 
the collapse the density increases 10 3 - 10 4 times. Since B oc p 2 / 3 , the magnetic 
field increases 10 2 - 10 27 times during the cloud formation (we are considering 
the magnetic field frozen in the plasma). For a pressure of ~ 5.5 x 10 4 dyn/cm 2 , 
the maximum magnetic field in the cloud, in order to have the magnetic field 
not dominating the pressure, is 

B 2 

— < 5.5 x 10 4 or B c < 8.3 x 10 2 G. 
4tt - 

In this way, the ambient magnetic field is then B a < (1.65 — 8.3) G. 

Our treatment can be understood in the light of the model of dos Santos et 
al. (1993a, b) that focuses on the mass loss from WR stars due to radiation 
pressure and Alfven waves. As radiation pressure line-driven models have 
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difficulties in explaining some observational characteristics of these stars (such 
as the disagreement between the observational mass loss rates and the maximum 
predicted mass loss by radiation pressure), they proposed that a flux of Alfven 
waves and radiation pressure act jointly. This fusion of an Alfven wave wind 
model (Jatenco-Pereira & Opher 1989a) and the radiation pressure CAK model 
resulted in good agreement with observations. In the present model we adopted 
the Alfven wave flux from the model of dos Santos ct al. (1993a, b), and showed 
that the heating due to the Alfven wave flux can cause a thermal instability 
which results in cold dense clumps coexisting with a hot diffuse gas. It is 
important to note also that near the star surface the magnetic force is more 
efficient than elsewhere (this is seen, for instance, in Figures 7 - 9 of dos Santos 
et al. 1993a). The thermal instability processes that we propose occurs just 
near the stellar surface. 

Another important aspect to discuss here is the competition between the 
instabilities acting in the wind. In principle, any thermal instability of the wind 
material would have to compete against the intrinsic line-driven instability of 
the flow. Then, consider that the cooling time, or thermal instability time, in 
the unperturbed medium, is 

tcooi ~ — 7— — l-6s (at Ph = 5.5 x 10 4 dyn/cm 2 ) , 
nA 

where T = Th, is the Boltzman constant, and A is the cooling rate 
dominated by radiative processes (Hb + H r {). Consider also that the dynamical 
time, or the line-driven instability time, at the base of the wind (up to one 
stellar radius), is 

R 

tdyn ~ — ~ 1.6 X 10 3 S , 

(v) 

where R is the distance from the ionization source (~ O.li?*) and (v) ~ 
5.04 x 10 7 cm/s (Owocki 1994; dos Santos et al. 1993a) is the wind velocity in 
this region. The above estimates show that the cooling time is much smaller 
than the dynamical time, as is necessary in order to have the thermal instability 
predominate in this region (see Krolik 1988; Mathews & Doane 1990). Moreover, 
the analysis of what instability is the most relevant for the thermodynamic 
of the wind is related to the region in which each instability operates. For 
instance, following Owocki (1994), we know that in the line-driven instability 
a minor part of the material is actually accelerated to high speed, for most of 
the mass the dominant effect is clumping. Diffuse radiation plays an important 
role in reducing the line-driven instability, especially near the wind base. The 
competition between these two instabilities may be important far in the wind, 
but at the base it is not. In this region thermal instability works alone. 

Thinking about the mass quantity involved in this condensation process, we 
can also note that if at a given pressure a thermal instability is indicated, namely 
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matter at a high temperature and a low temperature is permitted, little can be 
said about the fraction of the matter that is at the high or low temperature. For 
example, for the Crab nebula observations we have that almost all the matter 
is found to be contained in the filaments (Wilson 1971). These filaments are 
generally attributed to have formed by a thermal instability (e.g., Gouveia Dal 
Pino & Opher 1989). 

In spite of the presence of a gradient in the temperature, we did not consider 
the effect of thermal conduction in our model. We proceeded in this way because 
thermal conduction is extremely reduced in the gas, due to the presence of the 
magnetic field. This reduction occurs in the direction perpendicular to the field 
lines (Field 1965), and it is so important that even small magnetic fields are 
sufficient in eliminating this component of the thermal conduction (for some 
applications, see Begelman & McKee 1990 and McKee & Begelman 1990). 

In conclusion, our principal goal in this study was to explore whether or 
not a thermal instability, assisted by Alfvenic heatings, can play a role in the 
base of early-star winds. In spite of the fact that a number of simplifications 
were adopted in this first investigation, our results limit the pressure range for 
the existence of the two-phase equilibrium at the base of these winds. Our 
calculations indicate that a thermal instability may be a viable mechanism to 
form clumps in the winds of early- type stars. 
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Figure 1: The module of the heating/cooling processes, for logPy = 4.74 
(P H = 5.5 x 10 4 dyn/cm 2 ), F nl = F t = F A = 10 4 and t = 0.2 (as in Fig. 
2-5, where the thermal instability process occurs for Tl < T < Th with 
T L = 3 x 10 4 K and T H = 10 7 K). ' 
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Figure 2: The logT x logP diagram (H(T, n) w 0), with Fa = 10 2 : nonlinear 
and turbulent Alfvcnic heatings. 
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Figure 3: The logT x logP diagram {H(T, n) w 0), with F A = 5 x 
nonlinear and turbulent Alfvenic heatings. 
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Figure 4: The logT x logP diagram (H(T,n) w 0), with Fa = 10 4 : nonlinear 
and turbulent Alfvcnic heatings. 
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Figure 5: The logT x logP diagram (H(T,n) w 0), with Fa = 10 5 : nonlinear 
and turbulent Alfvcnic heatings. 
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